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This paper descibes a new method of calculating the mean duration and 
mean age of onset of a chronic disease from incidence and mortality rates. 
It is based on an ordinary differential equation resulting from a simple com- 
partment model. Applicability of the method is demonstrated in data about 
dementia in Germany. 
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1 Introduction 

This paper deals with analytical methods for calculating the mean age of onset and the 
mean disease duration in chronic diseases. For this calculation so-called IPM models 
are used. IPM stands for incidence, prevalence and mortality. First, the general IPM 
model is introduced. Then, formulas for the mean duration and the mean age of onset 
of the disease will be developed and applied to epidemiological data about dementia in 
Germany. 

Since we are interested in basic epidemiological parameters such as incidence, preva- 
lence and mortality with respect to a disease, it has proven to be helpful to look at 
so-called state models (compartmental models). The model that will be used here con- 
sists of the three states Normal, Disease, Death and the transitions between the states. 
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Normal means non-diseased with respect to the disease under consideration. The num- 
bers of persons in the Normal and Disease state are denoted as S (susceptible) and 
C (cases), respectively. The transition intensities (synonymously: rates) are called as 
shown in Figure [D z is the incidence rate, mo and mi are the mortality rates of the 
non-diseased and diseased persons, respectively. In general, the intensities depend on 
calendar time t, age a and sometimes also on the duration d of the disease. 




Figure 1: Simple model of a chronic disease with three states. Persons in the state 
Normal are healthy with respect to the considered disease. In the state Disease 
they suffer from the disease. The transition rates depend on the calendar time 
t, on the age a, and in case of the disease-specific mortality mi also on the 
disease's duration d. 



Henceforth, we assume that the rates do not depend on calenda r time and not on 
the duration of the disease. In this situation iMurrav and Lopea considered a two- 
dimensional system of ordinary differential equations (ODEs), which expresses the chang e 
in the numbers of healthy and sick patients aged a with the corresponding rates, (|l994l ): 



2 



— = - (i(a) + mo(a)) • 6 
da 

dC 

— — = 2(a) • S — mi{a) ■ C. 
da 



(1) 



Age plays here the role of temporal progression. In the paper (jBrinkd . l201ll ) this 
two-dimensional system is used to derive a scalar ODE of Riccati type, which relates the 
change in the prevalence at age a to the rates i, rriQ and mi: 



p)-{i-p - (mi - mo)). 



(2) 



In general, this Ricatti ODE has no closed analytical solution, we have to be contented 
with numerical solutions. However, there are special cases (e.g. mo = mi) in which the 
ODE is linear, and then there are explicit expressions for the solution. 

The benefits of such ODEs are obvious. For given age-specific incidence and mortality 
rates the age profile of the prevalence can be obtained by solving the Ricatti ODE. 
Typically this is called the forward or direct problem: we infer from the causes (the 
rates) to the effect (the prevalence). With smoothness constraints, the incidence and 
mortality rates determine the prevalence in an unambiguous way. 

In addition, for given prevalence and mortality rates the incidence can be dissolved. 
This is the inverse problem - we conclude from the effect to the cause. This allows, for 
example, cross-sectional studies being used for incidence estimates, for what otherwise 
lengthy f ollow-up studies would be necessary. An example of such an application is 



shown in (jBrinks , 



2011 
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2 Mean disease duration and median age of onset 



In this terminology, formulas for the mean age A of onset and the average duration D 
of a chronic disease can be derived. After the preparatory work in the introduction this 
is now quite simple: For the mean duration D the total person time spent in the disease 
state, i.e. the integral over C(a), has to be divided by the total number of new cases. 
The total number of new cases is the integral over the product of incidence and number 
of susceptibles. Hence, it holds 

JC{a)da fp{a) ■ N{a)da 

D = ^ = ^ , (3) 

/ i{a) ■ S{a)da f i{a) ■ (1 - p{a)) ■ N{a)da 



where N{a) := S{a) + C{a) is the total number of persons in the population with age a 
and UJ is an age that exceeds the age of the oldest member of the population. In inter- 
preting this expression as the first moment of a random variable D, the corresponding 
variance is 

J{C{a)-D) da J {p{a)-N{a)-D) da 

Var(D) = ^ = ^ . (4) 

/ i{a) ■ S{a)da J i{a) ■ (1 - p{a)) ■ N{a)da 



For the mean age A of onset, the new cases are weighted by the corresponding age. 

id LU 

J a ■ i{a) • S{a)da J a ■ i{a) ■ (1 — p{a)) ■ N{a)da 

A='-. ='-^ , (5) 

/ i{a) ■ S{a)da f i{a) ■ (1 - p{a)) ■ N{a)da 



Correspondingly, the variance of the age of onset is 

f {a-A) ■ i{a) • S{a)da / (a - A) ■ i{a) ■ (1 - p{a)) ■ N{a)da 

Var(.l) = ° = ° . (6) 

/ i{a) ■ S{a)da J i{a) ■ (1 - p{a)) ■ N{a)da 



Now we have these representations of D and A. Since the age-specific prevalence p 
is given by the ODE ([5]), D and A are completely determined by rates i,mo,mi and 
N. It is remarkable, that D and A depend on the shape N(a) of the age pyramid. One 
might expect that these numbers are characteristics that just depend on the disease, 
that they are disease inherent. However, easy considerations show that D and A have 
to be dependent on the age-distribution of the population. 
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The function N{a) is subject of demography. There are in the population models where 
the numbers of people in the age groups can b e represented analytically . The simplest 
example is the so-called stationary population, (jPreston and Coald . Il982l ) , which can for 
example be obtained from the following conditions: 



1. Mortality rates do not depend on calendar time. 

2. The birth rate does not depend on calendar time. 

3. At each age level, the net migration is 0. 

For a stationary population, there are explicit representations for N(a). However, 
real populations typically are non-stationary, which have to be managed differently. In 
Germany, there is an official population statistics from the Federal Statistical Office, 
which captures the age structure of the German population accurately. With N{a) 
given for every age a = 0, . . . , 100-1- from the official statistics, the integrals in Equations 
([3]) - ([6]) are replaced by sums. 



3 Application to dementia in Germany 



Our goal now is to apply the tools developed above to epide miological data about de- 
ment ia in Germany. The age-specific incidence is taken from (jZiegler and Doblhammer . 
2003 ). which reports data for males and females separately. The mortality m in the 
Germa n geri eral population is taken from the current life table of the Federal Statistical 
Office (|201ll ). Ref erence year is 20 10. The relative mort ality of the patients is set con- 
stant to R = 2.4, (jRait et al.l . l2010l ). as in (jBrinksl . l201ll ). In case the general mortality 
m and th e relative mo rtality R are given, the ODE ([2]) changes its type and becomes 
Abelian, (|Brinksl . \201Vi . Tab. 1): 



p)-{i-m- (l-[p-{R-l) + ir^^]. 



(7) 



Figure [2] shows the age-specific in cidence of dementia for males and females as reported 
in (jZiegler and Doblhammerl . l2009l ). For this work, the data have interpolated affine- 
linearly using the middle of the age classes as knots. 
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Figure 2: Age- and sex-specifi c incidence of dementia in Germany, 



( Ziegler and Doblhammeii . l2009l . Tab. 3) 
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The age-specific prevalence is derived by integrating the Abehan ODE ^ with initial 
condition p(60) = v i a the classical Runge-Kutta method (of fourth order), cf. e.g. 



(|Dahlquist and Biorckl . Il974l ) . All calculations are performed with the Software R (The 
R Foundation for Statistical Computing), version 2.12.0. The result is shown in Figure 

El 
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Figure 3: Age- and sex-specific prevalence of dementia after integration of the ODE d?]). 



Prevalence at age 60 years starts at 0, which is the initial condition. Until 70 years of 
age prevalence of men and women are almost the same. In the early 70ies the curves start 
to diverge, which is likely to be an effect of the different general mortality. Incidence rates 
in this age class are almost the same for men and women, however general mortality m 
for men is almost twice as high as for women. A higher general mortality mmaie > "^female 
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in Eq. ([7]) leads to a lower gradient ^^^^^'^ < ^ j™'"'" • It is striking that both prevalence 
curves have a maximum at age a* = 96 years. At this age it holds = 0, which implies 
i{a*) = p{a*) ■ mo(a*) ■ {R — 1). 

The age distribution of the new cases i{a)-{l—p{a))-N{a) for each age a = 50, ... , 100+ 
in 2010 is shown in Figure HI 
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Figure 4: Age distribution of the new cases (blue: males, red: females) in 2010 based on 
the age pyramid N{a) of Germany in the same year. 

The modus of the new cases is 80 and 85 years in men and women, respectively. 
Obviously, women are far more affected than men. The discontinuities stem from the 
discontinuous structure of the age pyramid. 
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The corresponding mean age A of onset and the mean duration D are presented in 
Table [D 

Table 1: Mean age of onset and duration of dementia in Germany in 2010. 

Mean age A Mean duration D 

(in years) (in years) 

Males 79^3 46 

Females 82.1 5.4 



4 Summary 

In this work a simple IPM model has been used to study the mean age of onset and mean 
duration of a chronic disease. These numbers depend on the incidence of the disease, the 
mortality of both diseased and non-diseased persons and also on the age pyramid of the 
population under consideration. The age-specific prevalence inherent in both numbers 
was ob tained as the solution of a new ODE, which has been derived in the predecessor 



paper (jBrinka . 120111 ). As a practical example, the calculations were applied to data 
about dementia in Germany. The mean age of onset of dementia is about 79 and 82 
years with mean duration 4.6 and 5.4 years in men and women, respectively. Due to 
the different life expectancies of men and women in Germany, the differences are not 
surprising. However, it is striking how big the difference in the nubmers of the cases 
of dementia are. Figure H] shows the differences in each of the age groups. In the year 
2010 a total of about 85000 men and 168000 women aged 65 years and above become 
dis eased. This is in the same order as the 78300 and 166000 males and females reported 



in (jZiegler and Doblhammerl . l2009l . Tab. 4) for the year 2007. The reasons for the about 
doubled number of women compared to men are on the one hand the higher incidence 
of dementia in females, which leads to a far higher prevalence (see Fig. [3|), and on the 
other hand the higher number of females aged 65-|-. For comparison, in Germany in 2010 
there are 7.2 and 9.6 million men and women aged 65 years and above, respectively. 

The methods described here are advantagous to predict characteristics of the persons 
with dementia. The age of onset might be important for comorbidities (such as diabetes) 
and associated late complications. The estimated number of diseased persons and disease 
durations are highly relevant for health services allocation planning. The methods allow 
predictions on regional levels, too. The age distributions of the states and communities 
in Germany differ quite substantially, which means that the functions N(a) in the states 
and communities are different. Hence the associated mean ages of onset and mean 
durations are different. 
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This study has some weaknesses. First, the method described needs the incidence 
and mortahty rates being independent from calendar time and disease duration. Both 
assumptions are hardly fulfilled. Hence, the calculations shown here are just an ap- 
proximation. Second, the incidences are based on claims data of the statutory health 
insurance (SHI) from 2002. Beside the fact that the data a re old compared to the ref- 



erence year 2010, SHI data tend to overestimate incidence, (j Abbas et al.l . I2OIII ). Third 



the ages equal to and above 100 are summarized in an age class 100-1- . The reason is the 
age pyramid of the Federal Statistical Office, which does not stratify ages beyond 100 
years. 

Another point is worth being mentioned: the examinations in this paper predict a 
peak in the prevalence of dementia in the second half of the ninth decade of life for men 
and women. After the peak, prevalence is decreasing again. If sufficient data is available 
for this age group, this hypothesis can be tested. 
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